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THROUGHPUT  AND  PACKET  ERROR  PROBABILITY  OF  CELLULAR 


FREQUENCY-HOPPED  SPREAD-SPECTRUM  RADIO  NETWORKS 


Jeffrey  W.  Gluck  and  Evaggelos  Geraniotis 


Department  of  Electrical  Engineering 
and  Systems  Research  Center 
University  of  Maryland 
College  Park,  MD  20742 


ABSTRACT 

We  characterize  multiple-access  interference  for  cellular  mobile  networks,  in  which  users 
are  assumed  to  be  Poisson-distributed  in  the  plane  and  to  employ  frequency-hopped  spread- 
spectrum  signaling  with  a  transmitter-oriented  assignment  of  frequency-hopping  patterns. 
Exact  expressions  for  the  bit  error  probabilities  are  derived  for  binary  coherently  demodulated 
systems  without  coding.  Approximations  for  the  packet  error  probability  are  derived  for 
coherent  and  noncoherent  systems  and  these  approximations  are  applied  when  forward-error- 
control  coding  is  employed.  In  all  cases,  the  effects  of  varying  interference  power  are  accu¬ 
rately  taken  into  account  according  to  some  propagation  law. 

Numerical  results  are  given  in  terms  of  bit  error  probability  for  the  exact  case  and 
throughput  for  the  approximate  analyses.  Comparisons  are  made  with  previously  derived 
bounds  and  it  is  shown  that  these  tend  to  be  very  pessimistic. 
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1  INTRODUCTION 

Research  in  the  areas  of  communication  networks  and  spread-spectrum  com¬ 
munications  has  failed  to  combine  the  two;  the  use  of  spread-spectrum  tech¬ 
niques  with  networking  techniques  in  order  to  provide  enhanced  multi-user 
capabilities  and  a  higher  degree  of  resistance  to  interference,  whether  hostile 
or  benign. 

Moreover  many  of  the  spread-spectrum  network  models  proposed  in  the 
past  have  lacked  precision  in  their  descriptions  of  the  effects  of  spread- 
spectrum  techniques  on  network  performance.  The  general  practice  in  pa¬ 
pers  like  [1]  and  [3]  has  been  to  model  the  effects  of  spread-spectrum  tech¬ 
niques  via  the  use  of  some  threshold  (a  number  of  users  beyond  which  com¬ 
munication  is  not  possible)  or  via  some  processing  gain.  Unfortunately,  this 
is  not  generally  considered  to  be  an  accurate  model,  and  the  results  obtained 
may  be  optimistic  or  pessimistic,  depending  upon  assumptions  made.  In  ad¬ 
dition,  these  analyses  do  not  accurately  account  for  data  modulation  and 
forward  error-control  coding. 

In  [1],  a  model  was  presented  for  a  frequency-hopped  (FH)  spread- 
spectrum  multiple  access  (SSMA)  digital  cellular  telephone  network  and  an 
expression  was  derived  for  throughput.  Notable  features  of  that  model  were 
the  characterization  of  user  mobility  in  terms  of  a  two-dimensional  Poisson 
process  and  the  treatment  of  transmitter  power  attenuation  with  distance. 
This  model  superceded  that  of  [18]. 

In  [2],  the  authors  described  and  analyzed  a  digital  cellular  mobile  net¬ 
work  using  some  of  the  techniques  of  [1],  While  adding  to  the  model  the  in¬ 
terference  due  to  users  in  other  cells,  the  authors  left  out  the  spatial  Poisson 
model  and  did  not  account  explicitly  for  power  attentuation  with  distance. 
Spatial  Poisson  modeling  was  also  used  in  [3],  in  which  the  authors  modeled 
spread-spectrum  multihop  networks.  In  [4],  a  more  precise  approach  was 
taken  to  the  spread-spectrum  multihop  network  problem,  while  using  some 
similar  network  features  in  the  analysis. 

In  this  paper,  the  modeling  techniques  described  above  are  combined 
with  a  precise  characterization  of  the  effects  of  frequency  hopping  for  the 
purpose  of  modeling  frequency- hopped  digital  cellular  networks.  This  anal¬ 
ysis  also  incorporates  different  data  modulation  types  (coherent  and  nonco¬ 
herent)  and  forward  error-control  coding.  Also,  unlike  several  of  the  previous 
analyses,  some  of  the  techniques  presented  here  allow  for  the  computation  of 
bit/symbol  or  packet  error  probability,  as  desired.  In  addition  to  modeling 
a  cell  of  a  FH/SSMA  cellular  mobile  radio  network,  this  also  corresponds 
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to  a  more  accurate  analysis  of  communications  in  the  presence  of  secondary 
multiple-access  interference  in  that  it  accounts  for  varying  powers  for  differ¬ 
ent  interferers  according  to  their  distances  from  the  receiver  (compare,  for 
example,  to  [11]  and  [19]). 

An  overview  of  the  contents  of  this  paper  is  as  follows.  In  Section  II, 
the  network  model  is  presented.  Section  III  contains  analyses  of  the  model 
for  coherent  and  non-coherent  modulation/demodulation  using  approxima¬ 
tion  techniques.  Section  IV  provides  an  exact  analysis  for  coherent  systems 
using  a  characteristic  function  method.  In  Section  V,  numerical  results  are 
presented  and  compared  with  previously  derived  results.  Section  6  contains 
the  concluding  remarks. 
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2  NETWORK  MODEL 


The  network  model  used  is  in  this  paper  similar  to  the  one  found  in  [1]. 
Users  transmit  packets  to  other  users  via  fixed  relay  (or  “base”)  stations. 
As  is  the  case  with  [1],  we  will  focus  on  the  user-to-station  channel.  Analysis 
of  the  station-to-user  channel  would  then  proceed  in  a  similar  fashion. 

Let  the  number  of  users  per  unit  area  and  hop  be  Poisson  distributed  with 
parameter  A.  Also,  let  M  =  (a:,  y )  represent  a  position  on  the  plane.  Denote 
the  origin  by  0  =  (0,0);  this  will  be  assumed,  without  loss  of  generality,  to 
be  the  position  of  the  fixed  station. 

Systems  like  this  generally  offer  voice  and  data  service.  In  the  case  of 
voice  service,  a  steady  stream  of  information  must  be  transmitted  and  delay 
can  not  be  tolerated;  however,  moderate  error  rates  can  be  tolerated.  In 
the  case  of  data  service,  traffic  is  bursty  and  large  delays  can  be  tolerated 
but  small  error  rates  are  necessary.  In  this  paper,  it  is  assumed  that  voice 
(telephone)  service  is  predominant. 

Slotted  ALOHA  is  taken  to  be  the  channel  access  protocol  and  the  sys¬ 
tem  is  assumed  to  utilize  (slow)  frequency  hopping  as  the  means  of  spreading 
the  spectrum.  Transmitter-oriented  assignment  of  hopping  patterns  is  also 
assumed.  Each  mobile  has  its  own  frequency- hopping  pattern  for  transmis¬ 
sion  and  reception;  the  base  station  has  a  list  of  all  of  the  patterns  and  can 
listen  to  several  of  them  simultaneously.  We  restrict  our  attention  to  the 
case  of  mobile-to-base  communication.  We  consider  the  situation  in  which  a 
mobile  station  has  already  established  communication  with  the  base  station 
at  the  center  of  the  cell  and'now  communicates  in  the  presence  of  secondary 
(or  multiple-access)  interference  from  other  mobile.  It  is  assumed  that  only 
the  other  terminals  in  the  cell  of  interest  may  interfere,  which  reflects  the 
use  of  other  frequency  bands  in  all  adjacent  cells,  and  that  frequencies  are 
only  reused  in  cells  distant  enough  so  that  the  effects  are  negligible. 

Synchronization  at  the  packet  level  is  assumed  feasible  for  all  users. 
Thus,  uncertainties  in  the  timing  between  different  users  can  be  small  rela¬ 
tive  to  packet  duration;  however,  since  they  might  not  be  small  relative  to 
the  dwell  time  of  the  frequency  hopper,  we  opt  for  a  general  asynchronous 
frequency  hopping  system. 


3 


3  APPROXIMATIONS  AND  BOUNDS  ON  PACKET 
ERROR  PROBABILITY  AND  THROUGH¬ 
PUT 


Our  approximation  methods  are  developed  below.  These  can  be  used  for 
coherently  and  noncoherently  modulated  systems  with  or  without  forward 
error-control  coding. 


3.1  Analysis  Based  on  Gaussian  Approximation  Techniques 


We  begin  with  noncoherent  signaling.  In  this  case,  the  system  is  assumed 
to  utilize  M- ary  frequency-shift  keying  (MFSK);  binary  FSK  is  a  special  case 
of  MFSK,  and  therefore  the  analysis  holds  for  this,  as  well.  Noncoherent 
demodulation  is  used  at  the  receiver. 

The  signal  transmitted  by  the  fc-th  mobile  user  is  of  the  form 

sk(t)=  n/2JS*(0  co${2ir[fc  +  bk(t)A  +  fk(t)]t+ek(t)},  (1) 


where  P0  is  the  transmitted  signal  power  and  is  assumed  to  be  the  same 
for  all  mobile  stations;  bk(t)  =  {2&j,^},  where  b 6  {1,2 ,...,M}  is  the 
k- th  sequence  of  M- ary  information  symbols  of  duration  T  each;  is  a 
pulse-shaping  waveform;  0k(t)  is  the  sum  of  phase  shifts  due  to  modulation 
and  frequency  hopping  generated  by  the  fc-th  transmitter;  fc  is  the  carrier 
frequency;  2 A  is  the  spacing  between  any  two  consecutive  tones  out  of  the 
M  different  tones;  and  fk(t)  is  the  fc-th  hopping  sequence. 

In  our  analysis,  we  consider  the  attenuation  of  signal  power  with  dis¬ 
tance.  Let  P(M)  denote  the  received  power  of  a  signal  originating  at  the 
mobile  station  located  at  the  point  M  of  the  plane  and  received  at  the  origin. 
We  will  assume  that  P(M)  is  given  by 


P(M)  = 


Por0°,  if  |0Af|  <  r0 
P0|OM|-“,  if  |0Af|  >  r0, 


(2) 


where  |0AfJ  denotes  the  Euclidean  distance  between  M  and  the  origin  0,  tq 
is  some  small  radius  within  which  P(M)  is  assumed  to  be  approximately 
constant,  Po  is  the  transmitted  power,  and  a  is  a  constant  that  reflects  the 
speed  of  attenuation.  r,o  is  necessary  because,  as  M  is  brought  arbitrary 
close  to  0,  |OM|-0'  approaches  infinity.  A  particular  signal  arriving  at  the 
receiver  (and  originating  from  Mk)  has  the  form 

rk(t)  =  yj2P(Mk)V(t)  cos  {2t r  [fc  +  bk(t)A  +  /*(*)]  t  +  £*(*)}  ,  (3) 
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where  \J2P(Mk)  represents  the  received  signal  amplitude  and  Ok(t)  is  the 
total  phase  shift.  Since  the  system  under  consideration  is  a  multiple- access 
one,  the  received  signal  has  the  form 

r(*)  =  r*(*  - T*)  +  n(*)>  (4) 

k 

where  each  rk(t)  is  given  by  (3),  rk  is  a  random  delay,  and  n(t)  is  zero-mean 
additive  white  Gaussian  noise  (AWGN)  with  two-sided  spectral  density 
Fig.  1  shows  the  receiver  for  the  i-th  signal  in  a  binary  FSK  system;  the 
M- ary  receiver  similar  in  form  but  with  M  branches.  The  ith  signal  is  the 
signal  of  interest.  r(t)  is  dehopped  and  the  dehopper  output  is  given  by 

"  V 

X  cos  [2-k  [fc  +  bk(t  -  rk)A]  +  <^(t)]}  +  nd(t).  (5) 


Here  /*(•)  and  /,(•)  represent  the  fc-th  and  ith  hopping  sequences,  respec¬ 
tively;  nd(t)  is  a  zero-mean  Gaussian  random  process  with  spectral  density 
6(u,v )  =  1,  if  u  —  v,  and  0,  otherwise;  and  <j>k(t)  represents  the 

total  phase  shift  of  the  signal  after  dehopping.  The  sufficient  statistics 
are  computed  as  shown.  Were  it  not  for  the  multiple- access  interference, 
the  computation  of  the  probability  of  error  would  be  straightforward,  as  for 
example  in  [10].  However,  this  difficulty  can  be  overcome  by  means  of  condi¬ 
tioning  the  error  probability  on  the  potential  number  of  interfering  users  and 
the  actual  number  of  interfering  users  causing  hits  on  the  desirable  received 
signal;  recall  that  a  hit  occurs  when  fk(t  -  rk)  =  fi(t)  (refer  to  [5]  and  [6]). 
To  evaluate  the  conditional  symbol-error  probability  just  described,  we  use 
the  Gaussian  approximation  technique,  as  developed  in  [5]  and  [6],  where 
the  Zc>m  and  ZSyin  were  approximated  by  zero-mean  Gaussian  random  vari¬ 
ables  with  the  same  second-order  moments.  Given  this,  the  symbol  error 
probability  can  be  derived  (see  [5]  or  [6])  and  is  given  by 


M-l 

ft,  =  £ 

m=l 


M-i  \  (-iy 

rn  I  m 


-f  1  6XP  [  2 (m  -f  l)a| 


m 


in  which  d\  is  the  variance  of  ZCiJn  and  ZJitn  and  is  given  by 


= 


(6) 

(7) 
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In  (7),  Eb  —  PqT  is  the  transmitted  signal  energy;  al  •  is  the  variance  of 
the  interference  with  the  t-th  signal  due  to  the  fc-th  signal;  and  g(|0M|)  = 
|0M|-°,  if  |0M|  >  ro,  and  Tq01,  if  |0M|  <  r0,  which  follows  from  (2). 

Regarding  the  effects  of  full  and  partial  hits,  their  independent  consider¬ 
ation  may  result  in  a  mathematically  intractable  expression,  if  we  attempt 
to  use  (6)  and  (7)  in  subsequent  expressions.  In  order  to  circumvent  this 
problem,  we  make  the  assumption  that  all  hits  are  full  hits.  Note  that  this 
causes  the  result  to  be  pessimistic,  since  we  assume  that  the  interference 
caused  by  any  hit  is  that  of  a  full  hit,  while  the  probability  of  a  hit  is  the 
sum  of  the  probabilities  of  full  and  partial  hits.  Under  this  assumption  (7) 
becomes 


a\  = 


|2  J^<7(|0M,|)log2  M  * 


+ 


X>(|0M*|),  (8) 


k*i 


in  which  aj  is  used  to  denote  the  variance  of  interference  under  the  assump¬ 
tion  that  any  hit  that  occurs  is  a  full  hit  (note  that  this  is  independent  of 
which  of  the  other  signals  is  interfering).  Expressions  for  aj  can  be  found  in 
[5],  in  which  it  is  also  shown  that,  for  random  FH  patterns,  o2  =  m^/M . 
For  a  rectangular  shaping  waveform  $(t),  m ^  =  1/3;  for  a  sine  shaping 
waveform,  =  (15  +  7r2)/127T2. 

At  this  point  we  use  a  three- term  exponential  approximation  to  Pc>j, 
the  symbol  correctness  probability,  an  exponential  approximation  is  required 
to  facilitate  the  remainder  of  this  analysis.  Such  an  approximation  can  be 
made  using  techniques  like  Prony’s  method  (see  [8]),  or  curve-fitting.  A 
least-squares  curve  fit  was  used  to  obtain  the  results  shown  below;  as  a 
result,  we  have 


Pc,i  —  1  —  Pe,i(x) 

3 

«  CuelvX,  (9) 

v-\ 

where  we  define  x  to  be  a\.  In  using  curve-fitting  for  (9),  we  considered 
the  range  of  x  to  be  the  interval  [0,2];  this  was  chosen  on  the  basis  of  the 
various  parameters  characterizing  the  statistics  of  user  population  and  the 
varying  distances  from  the  base  station.  If  we  substitute  for  x  in  (9)  and 
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substitute  (8)  for  o\,  we  get 


Pc,i  ~  Yj  C"  eXP 


l/=l 


lv  I  |2^/7(|°M,|)log2M 


+ 


Let  us  now  write  (10)  in  the  form 

3 


pc,i  =  Yc‘'n(v’f)eSv{f)^k*i9{[0Mk[)  ■ 


(10) 

(11) 


i/=i 


In  (11),  we  have  introduced  the  notation 


exp 


^■jfQ9{r)\og2  M 


and 


K{f) 


7i/<7/ 

9(f) 


in  order  to  simplify  this  and  subsequent  expressions.  Also,  we  have  taken 
f  =  |0M<|,  the  distance  of  the  “terminal  of  interest.” 

Let  us  reiterate  that  the  expression  given  in  (11)  represents  the  con¬ 
ditional  probability  of  symbol  correctness,  which  is  conditioned  on  f  (the 
number  of  potential  interferers  present),  and  on  the  number  of  hits  that  ac¬ 
tually  occur  given  the  above  number  of  potential  interferers.  However,  since 
our  ultimate  goal  is  to  determine  packet  throughput,  we  require  a  packet 
error  correctness  probability. 

The  development  of  the  probability  expression  that  corresponds  to  (11) 
above  for  the  coherent  case  is  similar  to  the  one  carried  out  above  for  the 
noncoherent  case.  In  this  case,  the  system  is  assumed  to  utilize  binary  phase- 
shift  keying  (BPSK)  with  coherent  demodulation.  The  signal  transmitted 
by  the  fcth  mobile  station  is  of  the  form 


sk(t)  =  y/2Pobk(t)'$(t)  cos [271- [/c  -f  fk(t)]t  +  0k(t)) ,  (12) 

in  which  all  quantities  are  as  above,  except  that  bn’  6  {— 1,+1).  We  now 
have  information  bits  rather  than  M- ary  symbols.  From  (12)  we  obtain  the 
following  expression  for  the  fc-th  received  symbol 

rk(t)  =  ^2 P(Mk)bk(t)V(t)  cos  [2tt  [fe  +  /*(«)] t  +  **(<)]  •  (13) 
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The  received  signal,  given  by  (4),  goes  through  the  same  dehopper  and 
the  dehopper  output  is  given  by 

rd(i)  =  £  | 

cos  [2irfct  +  <f>k(i)])  +  MO-  ( 14) 


iP(Mk) 


6[fk(t  -  Tk),fi(t))bk(t  -  Tfc)®(t  -  n) 


This  is  fed  into  a  coherent  demodulator  and  the  sufficient  statistic  for  the 
reception  of  the  Zth  bit  of  the  tth  transmitted  signal  is  given  by 

(f+l)T 

Z,  =  J  rd(t)^/(t)  cos(27r fct)dt.  (15) 

IT 


Assuming  that  the  receiver  is  synchronized  to  the  ith  signal,  Z,  can  be 
rewritten  as 


Zi  = 


T  tj0 


+  £ 


P(Mk) 


T  h,i  4-  n, 


(16) 


where  Ik,i  represents  the  normalized  interference  caused  by  the  fcth  signal 
to  the  reception  of  the  ith  signal  and  n  is  a  zero-mean  Gaussian  random 
variable  with  variance  N0T/16. 

The  decision  process  is  concerned  with  determining  whether  the  current 
bit  is  +1  or  —1.  The  test  is,  therefore,  to  determine  whether  Z;  >  0  or 
Zi  <  0.  This  allows  some  further  manipulation  of  Zi.  Zi,  the  modified  test 
statistic  after  we  divide  by  P0,  is  of  the  form 


Zi  =  V^CIOM.U&W  +  ^2  \] g{\0 Mk\)h,i  +  n,  (IT) 

k*i 


where  n  is  a  Gaussian  random  variable  with  variance  No/2Eb2Eb,  where 
Eb  =  PqT  is  the  transmitted  bit  energy.  We  utilize  the  Gaussian  approx¬ 
imation  technique,  namely  we  approximate  Z,  by  a  zero-mean  Gaussian 
random  variable  with  variance  ,  which  denotes  the  variance  of  Z,.  This 
is  the  coherent  analog  to  the  Gaussian  approximation  found  in  [5]  and  [6]. 
Given  this  approximation,  the  conditional  symbol  (bit)  error  probability  is 


Pr,i 


-  erfc 

2 


yjpg 


(18) 
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in  which  erfc(-)  is  the  usual  complementary  error  function.  £From  our  earlier 
expressions,  it  can  be  shown  that 


a2,  =  \  If  +  S  9(\0Mk\)<rli,  (19) 


where  a\  i  is  the  variance  of  7^,.  If  we  once  again  assume  that  all  hits  that 
occur  are  full  hits,  we  find  that 


1 

2 


'TV 

No. 


k^i 


(20) 


where  aj  is  the  variance  of  the  interference  terms  under  the  assumption  that 
any  hit  that  occurs  is  a  full  hit.  It  can  be  shown,  using  the  results  of  [5], 
that  <Ty  =  m^,,  where  m ^  is  as  'discussed  above. 

At  this  point  we  use  a  three-term  exponential  approximation  to  Pc  t  to 
obtain 


Pci  =  1-^  =  1-- 


£  c- 


(21) 


i/=i 


where  we  define  x  to  be  2ct^  /<7(|0M;|).  Again  we  use  [0,2]  as  the  range  of 


x.  If  we  substitute  for  x  in  (21)  and  substitute  (20)  for  cr|  ,  we  get 


Pc,i  ~  Gv  exp  (7^ 
2  m^, 


U=1 


(KIOMfll)'1 


+  S(|0M,l)  ^ 

H(v,r)  and  6u(f)  are  given  by 

H{v,f)  =  exp  I7*  (ff(f)^) 


(22) 


-1 


and 


6y{f) 


27 

9(r) 
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This  completes  the  development  of  the  conditional  coherent  bit  correctness 
probability. 

As  mentioned  earlier,  in  order  to  compute  packet  throughput  we  need  a 
probability  of  packet  correctness  or  error.  We  would  like  to  get  an  expected 
packet  error  probability  with  respect  to  r  and  the  number  of  users.  We  do 
this  by  taking  a  known  error  probability  expression  and  substituting  it  in 
our  result.  Such  probability  expressions  can  be  found  in  or  adapted  from 
expressions  found  in  many  sources  (e.g.  [9]).  Expectations  are  then  taken 
with  respect  to  f,  the  number  of  users,  and  the  number  of  hits.  This  method 
is  demonstrated  below  for  the  cases  of  uncoded  packets  and  RS(n,fc)-coded 
(encoded  with  a  Reed-Solomon  code  of  length  n  and  rate  k/n )  packets;  the 
development  is  the  same  for  other  coding  schemes. 

For  a  general  uncoded  packet  of  length  L  we  have 

Pe,vkt  =  1  -PCL. 

Since  Pc(r,  Jr)  in  (11)  is  dependent  on  r,  on  the  number  of  potential  inter¬ 
feres,  and  on  the  number  of  hits  given  K ,  we  must  take  expectations  over 
all  of  these.  Thus, 

P',?kt  =  1  -EfEKEh{PcL(r,K)}.  (23) 

The  inner  expectation  in  (23)  requires  some  explanation.  The  probability 
of  having  k  hits,  when  there  are  in  addition  to  the  user  i,  K  other  users 

P£(  1  —  Ph)K~k,  where  Ph  is  the  probability  of  a 

full  hit.  As  discussed  in  [7],  for  stationary  Markov  random  hopping  patterns 
Ph  is  given  by 

where  q  is  the  number  of  different  frequencies  and  N3  the  number  of  symbols 
per  hop.  As  discussed  in  [7],  for  moderately  large  values  of  q,  (24)  either 
bounds  or  approximates  the  hit  probabilities  for  memoryless  random  hop¬ 
ping  patterns  and  Reed-Solomon  periodic  hopping  patterns.  Thus  (24)  is  a 
good  general  expression  for  most  purposes. 

Suppose  now  that  the  summation  in  the  exponent  in  (11)  can  be  decom¬ 
posed  into  collections  of  interferers  having  equal  g(jOM|).  We  would  then 
have 

]£  S(|0A**I)  =  2C  k>9P  (25) 

k^i  j= 1 


present,  is  given  by 
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where  kj  is  the  number  of  interferes  with  5(jOM|)  =  gj. 

Let  us  consider  some  planar  region  A  and  a  partition  of  A  given  by 
such  that  for  all  mobile  users  located  in  Aj,  5(|0M|)  =  gj.  Let 
N(Aj)  denote  the  number  of  such  mobile  users.  We  now  have  a  total  number 
of  users  with  a  given  attenuation  and  a  probability  that,  out  of  a  given 
population  of  K  potential  interferes,  there  are  k  hits.  Then  one  may  take 
Kj  =  N(Aj)  to  be  the  population  in  that  probability  expression. 
iFrom  (11)  we  get  by  expanding  the  trinomial 

phf.K)  = 

m=0  t=0  \  J  \  / 

x  [C3H( 3,  exp  +  iS2(f) 


+(L  -  m  -  t)63(f)]^p(|0M^|) 

fc/i 


(26) 


Considering  this,  the  discussion  oiEh{-},  and  the  fact  that  C,(f)  =  CiH(i,r), 
we  remit  (23)  as 


P  e,pkt  — 


L  m 

i-EE 

m=0  i=0 


L 

m 


m 

i 


xE?EKEh  {(51m(f)Ci(f)C3L~m_‘(0 

x  exp([m6i(f)  +  i62(f)  +  {L  -  m  -  1)63(7=)] 


X  Yj9ik: 
i=i 

L  m  / 

i-EE 


i=0  t=0  \ 


L 

m 


m 

i 


x£,{cr(f)C5(f)C3£'-"'-'(f) 


xEk 


lasts) 


Pkh’(l-Phf>-h’ 


X  exp  ((m^i(f)  +  iS2(f )  +  (L  -  m  -  z)63(r)] 

v  / 

L 


L  m 

■-EE 

m= 0  i=0 


m 


m 

i 
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xEr{c?(f)Ci(f)Ctm-\?) 

U=i 

+1-Ph]*'}} 


L 

m 


m 

i 


L  m  / 

i-EE 

m=0  i=0  \ 

xPr-{cr(0Q(^3_m"‘(r) 

xEk 


|  Y1  Ki ln  [Aexp([m^i(f)  +  iS2(f) 
+(I  -  m  -  i)63{f)]gj)  +  1  -  Pft]}}  . 


(27) 


Ek{ •}  is  taken  by  following  a  technique  used  in  [1].  For  any  partition  of 
A,  say  {Aj}”-!,  and  constants  {a,}"=1, 


E  <  exp 


Li=1 


=  exp 


L  i=i 


(28) 


where  S(Aj )  denotes  the  area  of  Aj.  Because  Kj  =  N(Aj) ,  (27)  and  (28) 
can  be  combined  to  yield 


L  m  / 

T'.*t  =  i-££( 

m=0  »=0  \ 


L 

m 


m 

i 


x£,{(;r(f)<5;(f)C3L-’"-i(f) 

n 

x  exp  A  ^2  S(Aj )  (P/,  [exp  ([ m6i(f )  +  t<52(f) 
.  i=1 

+(L-m-  i)S3{f)]gj)-  1])]}. 


(29) 


Generally,  <7(|0Af|)  can  be  considered  to  be  a  limit  of  step  functions,  which 
allows  Pe<pkt  to  be  written  as 


L  m 

Pe,pkt  —  1  _  Xy 

m=0  i=0 


L 

m 


m 

i 
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X  exp  j  ( Ph.  [exp([m^i(f)  -f  iS2(f ) 

+(L-m-i)63(r)]ff(lOM\))-l])]}  (30) 

where  fA(-)d2M  denotes  the  integral  over  the  surface  A.  This  completes  the 
analysis  because  Ef{-}  is  simply  an  expectation  over  the  entire  area  A  under 
consideration. 

We  now  consider  RS(n,A:)-coded  packets  with  error-only  correction  de¬ 
coding.  The  applicable  error  probability  expression  is  given  by 

P',rs  =  i  -  £  ( " )  cr'(>  -  ft)',  (31) 

where  t  —  |_n  —  k( 2j  is  the  error  correction  capability  of  the  code.  This 
can  also  be  expressed  as 


^=si0t)(-ir^<+"  (32) 

With  respect  to  the  required  expectations,  this  can  be  treated  exactly  like 
the  uncoded  finally  packet  error  probability  expression.  Therefore,  we  obtain 


Pe,RS 


t  t  n—l+m  i 

‘-EE  £  £ 

/~0  m= 0  »=0  ;'s 0 


x  /  *  j  (-1)"*  {q(f)<;i(f)c3”-'+"-i-i(f) 

x  exp  |a  (Ph  [exp  +  jS2(f) 

+(«  -  l  +  m-i-  j)63(f)]^(|0M|)) 
-l])d2M]} 


(33) 


It  should  be  noted  that  in  the  case  of  coherent  BPSK,  the  three-term  expo¬ 
nential  approximation  must  approximate  P™-,  rather  than  PCi,-  when  consid¬ 
ering  Reed-Solomon  codes.  The  m  noted  here  refers  to  the  number  of  bits 
in  BPSK  per  Reed-Solomon  symbol;  for  example,  if  RS(32,16)  code  is  used, 
m  =  5.  It  should  also  be  noted,  when  discussing  Reed-Solomon  codes,  that 
the  signal-to-noise  ratio  is  normalized  by  the  code  rate. 

Once  one  has  obtained  an  expression  for  packet  error  probability,  the 
expected  throughput  can  be  obtained  by  subtracting  it  from  one  and  mul¬ 
tiplying  by  the  offered  traffic  load,  which  is  given  by  Gtoi  —  Att5(v4), 
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S(  A)  being  the  area  of  the  region  of  concern.  In  the  cases  presented  here,  a 
bandwidth  expansion  factor  must  also  be  introduced. 


S  =  Gtot(l-Pe)B. 

(34) 

B,  the  bandwidth  expansion  factor,  is  given  by 

B  =  - Rc , 

<7 

(35) 

where  q  is  the  number  of  hop  frequencies  and  Rc  the  code  rate  Rc  is  1 ,  for 
uncoded  packets,  and  k/n,  for  RS(n,k)  coded  packets. 


3.2  Bounds 

Error  probabilities  have  been  upper  bounded  giving  lower  bounds  on  system 
throughput.  These  will  be  presented  here  for  the  purpose  of  comparison. 
The  probability  of  a  symbol  error  at  the  receiver  p,(K),  when  K  users 
transmit,  can  be  upper-bounded  as 

p,(A0  <  1  -  (1  -  Po)(l  ~  (36) 

where  P0  is  the  error  probability  for  a  single  user  system  in  the  presence 
of  AWGN  with  or  without  fading,  and  Ph  is  the  probability  of  a  hit  from 
a  single  other  user;  expressions  for  Ph  can  be  found  in  [7].  For  random 
frequency-hopping  patterns,  we  can  use  Ph  =  (1  +  log2  M/iVfc)  where  Nb 
denotes  the  number  of  bits  per  dwell  time  and  q  is  the  number  of  frequencies. 

Let  us  discuss  now  the  calculation  of  throughput  when  the  total  number 
of  users  (packets)  transmitting  in  a  slot  (a  packet  slot,  which  is  L  symbols 
long  in  the  uncoded  case)  is  K.  Let  s(K)  and  Pe(K)  denote  the  average 
number  of  packets  received  successfully  and  the  average  probability  of  failure 
(erroneous  reception)  of  a  typical  packet,  respectively,  when  K  packets  are 
transmitted.  This  gives  use  s(K )  —  K  —  KPe(K). 

Let  K  be  Poisson-distributed  on  an  area  A  with  average  intensity  A 
(packets  per  unit  area).  Let  se(G )  and  Pe(G)  denote  the  averages  of  the 
quantities  s(K)  and  Pe(K)  defined  in  the  previous  paragraph  with  respect 
to  the  Poisson  distribution  p(  K)  =  e~GGK /K\,  where  E{K)  =  A  A  —  G 
(the  offered  load).  Then  Pe(G)  =  1  -  se(G)/G  =  E{KPe(K)}/G  and 
the  normalized  throughput  is  given  by  q(G)  =  G/q  ^1  -  Pe(G')j. 

We  can  actually  show  that  Pe(G )  =  E  {Pe(K  +  1)}-  Since  the  packet 
error  probability  Pe(K )  is,  through  the  symbol  error  probability  P,{K ),  a 
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function  of  K  —  1  (the  number  of  interfering  packets),  PC{G)  is  a  function  of 
K  and  not  of  K  —  1.  This  implies  that  in  the  calculation  of  the  throughput, 
instead  of  working  with  the  total  number  of  transmitted  packets,  we  can 
equivalently  work  with  the  number  of  interfering  packets  and  treat  that  as 
being  Poisson-distributed  with  the  same  parameter  G.  This  argument  also 
justifies  the  assumption  that  the  number  of  interfering  users  in  Poisson- 
distributed  with  parameter  G. 

We  can  now  write 

PC(G)  =  1  -Pe(G)  =  E  {[1  —  pa(K  +  1)]}L 

<  je|[(i-p0)(i-p,)^]£'} 

=  (37) 

where  L  is  the  number  of  Af-ary  symbols  or  bits  transmitted  in  a  packet 
L  =  1  can  be  used  to  bound  the  bit  error  probability.  The  throughput  can 
be  then  lower-bounded  as 

r j(G)  >  ^(1  -  P0)VG[1-<1“Pfc>i].  (38) 

If  the  system  employs  error-control  coding,  we  can  follow  the  analysis 
in  [11].  The  analysis  will  be  shown  for  a  Reed-Solomon  (RS)  errors-only 
correction  scheme;  other  schemes,  such  as  RS  with  side  information,  have 
been  analzyed  similarly. 

The  probability  of  a  symbol  error  ps{K)  is  upper-bounded  in  a  way  sim¬ 
ilar  to  that  in  REF(36),  where  Ph  is  the  probability  of  a  hit  for  an  RS  sym¬ 
bol,  consisting  of  mc  m- ary  symbols;  thus,  we  must  use  mc  log2  M  instead 
of  log2  M  in  the  expression  that  provides  the  probability  of  a  hit.  When 
the  RS(n,k)  code  errors-only  correction  is  employed  and  one  codeword  per 
packet  is  transmitted,  the  probability  of  correct  reception  of  a  codeword  is 
given  by 

Pc(K)  >  £  (  "  )  )]J  [1  -  pMO)n-j,  (39) 

where  t  =  [n  —  fc/2l  is  the  error-correction  capability  of  the  code  and 
P»(K)  <  1  —  (1  —  P0)  -1.  Finally,  we  can  get  an  expression  for  the  aver¬ 
age  probability  of  correct  packet  reception  PC(G )  that  allows  us  to  write 
the  following  expression  for  the  normalized  throughput  (per  frequency  slot) 
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V  (G): 

v(G)  >  “£of^E(-r)[i-(i-A,r(i-p^]J 

x  [(l-Po)m‘(l-Ph)K]n~3 

=  ~tt(")(i)  (-i)'(i  -  p0)"‘.(«+«w).-o[*-a-n)-+(4fl) 
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4  ACCURATE  EVALUATION  OF  BIT  ERROR 
PROBABILITY 

If  coherent  modulation  is  employed,  the  resulting  system  can  be  analyzed 
by  means  of  a  characteristic  function  method,  as  will  be  demonstrated  below. 

4.1  Analysis  Based  on  Integrating  the  Characteristic  Func¬ 
tion  of  the  Interference 

In  this  case,  the  system  will  be  assumed  to  utilize  binary  phase-shift  keying 
(BPSK)  with  coherent  demodulation.  The  signal  transmitted  by  the  fcth 
mobile  station  is  of  the  form  given  in  (12).  We  assume  here  that  Rayleigh 
fading  may  or  may  not  be  present.  We  modify  (13)  to  get 

nt(0  =  Fk^2P(Mk)bk(t)^(t)  cos[27r[/c  +  Jk(t)]t  +  0*(<)],  (41) 

where  Fk  is  a  Rayleigh  fading  factor,  if  fading  is  assumed  to  be  present  and 
all  other  quantities  are  as  previously  described. 

The  dehopper  output,  rd(t)  which  is  used  in  (15)  is  now  given  by 


X  cos  [27t  fct  -f  <pk(t)}}  +  nd(t).  (42) 

We  proceed  until  we  again  have  a  modified  test  statistic  Zi  of  the  form  given 
in  (17).  Our  goal  is  to  derive  expressions  for  the  probability  of  error  and 
the  expected  throughput  r/.  Since  a  coherent  receiver  is  being  used  here,  the 
characteristic  function  method,  described  in  [13]  and  simplified  to  the  case 
where  the  signature  sequences  are  identically  1,  will  be  used  to  derive  the 
desired  probability  expression. 

The  computation  of  the  error  probability  involves  the  characteristic  func¬ 
tion  of  the  noise  n,  which,  in  our  case,  is  $^(u)  =  exp  (— u2/ 2  •  No/2Ebto), 
and  the  characteristic  function  of  the  interference  process,  which  is  the  mid¬ 
dle  term  in  (17).  This  will  be  derived  here  for  the  synchronous  case  and  for 
the  asynchronous  case,  both  with  and  without  Rayleigh  fading. 

In  the  synchronous  frequency-hopping  case,  we  have  Tk  =  0  for 
all  k.  The  fc-th  signal  interferes  with  the  t-th  signal  if  and  only  if  the  two 
signals  are  exactly  superimposed  upon  each  other;  this  suggest  that  user  k 
interferes  with  user  i  if  there  is  a  full  hit,  when  two  signals  are  sent  in  the 


rd(t)  =  ]£ 


>P(Mk) 


b[fk{t  -  rk),  fi(t)}bk(t  -  rjt)$(t  -  rk) 
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same  frequency  “bin”  at  the  same  time.  Thus,  if  q  denotes  the  number  of 
different  frequency  bins,  the  probability  of  a  hit  is  1  jq. 

A  result  used  in  [1]  can  now  be  adapted  to  give  an  expression  for  $/(u), 
the  characteristic  function  of  the  interference  process.  The  result  used  in  [1] 
is  that,  if  the  interference  process  can  be  written  as  l  —  Ylk  Xkh(Mk)  for  a 
random  process  {X*;}  with  pdf  px(-)  and,  if  the  number  of  users  is  Poisson 
distributed  on  the  area  A  with  parameter  (mean  intensity)  A,  then 


<h/(u)  =  exp  A  J  ^  J  px(x)e,uxh^M^dx  -  1^ 


d2M 


(43) 


in  which  fA  f(M)  d?M  denotes  a  surface  integral  and  A  represents  the  area 
of  the  cell  (area  of  concern).  This,  however,  is  not  the  exact  form  that  must 
be  used  here,  since  the  interference  process  in  our  problem  has  the  form  of 
the  middle  term  of  (17)  and  the  probability  of  hits  must  also  be  taken  into 
account.  As  a  result,  the  characteristic  function  of  the  interference  process 
takes  the  form 


$/(u)  =  exp  ^  -  1  d2M  j,  (44) 


where 

1  f2* 

4>(u)  =  —  /  exp(tucos0)d0  =  Jo{u).  (45) 

2tt  Jo 

In  (45),  the  cos 6  term  was  obtained  by  considering  (43)  and  (17),  letting 
Tk  =  0  and  assuming  uniformly  distributed  random  phase  processes  (6  is 
uniform  on  the  interval  [0, 27t));  Jo(-)  is  the  Bessel  function  of  the  first  kind 
and  zeroth  order. 

If  we  consider  our  region  of  concern  to  be  a  disk  of  radius  R ,  then  the 
integral  in  (44)  can  be  easily  computed.  Let  us  denote  by  r  the  radial 
distance  \OM\  and  consider  the  annulus  (r,  r  +  dr).  Then  we  can  write  that 
Sa  [*Jo  (uv^flOAfl))  -  lj  d2 M  =  /(f  2nr  [j0  (uv/fi,(r))  ~  l]  dr,  where  g(-) 
is  as  in  (17).  This  expression  can  be  evaluated  using  formulas  found  in  [14] 
to  yield  series  solution. 

Let  us  consider  the  same  model  with  Rayleigh  fading.  This  is  accounted 
for  by  means  of  Fk  in  the  expressions  above.  The  Fk  are  assumed  to  be 
independent  and  identically  distributed  with  a  Rayleigh  parameter  b 2.  As  a 
result  of  the  addition  of  fading,  (45)  becomes 


$(u) 


exp(tuy  cos  6)d6dy 
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(46) 


) 


In  (46),  Pf{v)  =  i/2/62exp  (— y2/262)  is  the  Rayleigh  density  function 
with  parameter  62.  Since  the  product  of  a  Rayleigh  random  variable  and 
the  sine  or  cosine  of  a  uniform  on  [0,27r]  random  variable  yields  a  Gaussian 
random  variable,  G  —  FcosO  has  pdf  Pg(j?)  =  l/y/2nbexp  (— rj2/262) 
and  the  right-hand  equality  is  derived  by  performing  the  integration  with 
respect  to  q. 

If  the  circular  region  model  mentioned  above  is  used,  (44)  becomes 


where  B{u)  —  b2u2  / 2.  A  closed-form  result  can  be  obtained  using  a  formula 
from  [16],  if  a  is  known.  If  the  typical  value  a  =  3.5  is  used,  (47)  becomes 

\7ne-B^x 


+ 


B{u) 


+  £(-I)  PwF . 


(48) 


We  now  focus  on  the  asynchronous  frequency-hopping  case.  It  will 
be  assumed  that  the  in  (42)  are  independent  and  identically  uniformly 
distributed  on  the  interval  [0,  T/,).  As  a  result  of  these  delays,  there  are  three 
possible  interference  scenarios  for  a  signal  transmitted  in  a  given  frequency 
bin:  no  hit,  partial  hit,  or  full  hit.  These  must  all  be  accounted  for  in  the 
computation  of  the  characteristic  function  of  the  interference  process. 

Since  the  interference  term  in  (17)  is  a  sum  of  independent  random 
variables,  the  overall  characteristic  function  can  be  seen  as  a  product  of  the 
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individual  characteristic  functions,  that  is, 


In  (49),  E  denotes  expectation  with  respect  to  the  number  of  interfering 
users,  Pj  is  the  probability  of  a  full  hit,  Pv  is  the  probability  of  a  partial 
hit,  1  —  P  =  1  -  Pj  -  Pp  is  the  probability  of  no  hit  occurring,  and  $^k,l\u) 

and  $pk,'\u)  are  the  characteristic  functions  of  7fc,j,  given  that  a  full  hit  or 
a  partial  hit  has  occurred.  Note  that  for  the  case  of  no  hit,  7^,-  =  0,  which 
results  in  a  characteristic  function  equal  to  1.  The  validity  of  the  second 
equality  in  (49)  was  proven  in  [13]. 

The  quantity  inside  the  expectation  in  (49)  can  be  written  in  the  form 
exp  |]Cfc/t'ln  [l  —  P  +  Pj^(k<,){u)  -f  Pp$pA:,^(u)J  j.  Let  us  consider  a  step 
function  g(\OM\)  =  5Zj=i  9j  where  is  a  partition  of  the 

region  A;  7b(A7)  is  1,  if  M  £  B ,  and  0  otherwise;  and  g(\OM\ )  =  g:  for 
all  M  e  Aj.  If  we  denote  by  N(Aj )  the  number  of  users  in  the  area  Aj,  we 
obtain  that 


(50) 


in  which  $/(uv/ffj)  and  ^(uy/^J)  denote  the  characteristic  functions  of  the 
interference  due  to  any  particular  user  located  in  Aj  for  the  cases  of  full  hits 
and  partial  hits,  respectively.  Then  the  expectation  can  be  written  as 


$/(u)  =  exp  |  XPj  ^2  S(Aj)  [^/(u-y/ffj)  -  1] 

;'= l 

In  the  general  case,  p(|OM|)  is  a  limit  of  step  functions,  implying  that 
$ /(u^/gdOMD)  and  ^(tt^dOMl))  are  such  limits,  too,  and  the  overall 
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characteristic  function  becomes 


$/(u)  = 


exp  |  A Pj-J  (uy/g(\OM\)^  -  l]  d2M 

+  XPp  jf  [$„  (u^/fldOAfl))  -  l]  d?M }  . 


(52) 


Expressions  for  Pj  and  Pp  were  given  in  [13]  for  first-order  Markov  ran¬ 
dom  hopping  patterns,  for  memoryless  random  hopping  patterns,  and  for 
Reed-Solomon  periodic  hopping  patterns.  Following  [13]  for  the  case  in 
which  the  signature  sequence  is  identically  1,  we  obtain  that 


*P(tO 


~  jT  Jo[^[R^r)  +  R^T)}^ 

2 7  /:  h  (ia*(r)) + ( 


dr 


dr , 


(53) 


in  the  absence  of  fading,  and 

* / iu)  =  f  Jq  expj-^- [R4l(r)  +  Rip(r)]2 


=  2?  Jo  {exp 

+  exp 

-b^-R%(r) 

dr 


(54) 


in  the  presence  of  Rayleigh  fading.  In  (53)  and  (54),  the  functions  R y,  and 
R ^  are  the  partial  autocorrelation  functions  of  the  shaping  waveform  ip(t) 
and  are  defined  by  R^(s)  =  jj  -  s)dt  and  R^(s)  =  R^(T  —  s). 

The  relationship  in  [13],  which  expresses  the  bit  error  probability  of  the 
receiver  as  an  integral  involving  the  characteristic  functions  of  the  noise  and 
the  interference  process  (see  [6]),  can  be  modified  to  depend  on  the  position 
of  the  “transmitter  of  concern”  in  the  following  way 

n(|0M|)  =  Q  ^Eb'o/No'j 

1  sin[j(|0Mj)u]  ,  .  r,  _  .  ,  . 

H - /  - -  $r»(u)  [1  -  $/(u)]  du,  (55) 

7T  Jo  u 

where  |0M|  represents  the  distance  between  the  receiving  (base)  station  and 
the  “transmitter  of  concern”  (the  ith  user  in  the  analysis  above),  Q(x )  is  the 
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complementary  error  function  given  by  Q(x )  =  exp  (— y2/ 2)  dy, 

and  $n(u)  was  defined  earlier.  $/(u)  is  given  by  (44)  for  synchronous  sys¬ 
tems,  by  (45)  for  no  fading,  and  (46)  for  Rayleigh  fading.  For  asynchronous 
systems,  it  is  given  by  (52)  and  (53)  for  no  fading  and  by  (54)  for  Rayleigh 
fading. 

The  result  in  (55)  is  useful  because  it  gives  us  a  measure  of  the  correctness 
of  each  individual  bit.  Unfortunately,  to  evaluate  the  throughput  one  needs 
the  packet  error  probability  Pe(|0M|),  which  can  not  be  expressed  as  an 
explicit  function  of  Pb(\OM\).  The  error  probability  of  each  one  of  the  L 
bits  in  a  packet  depends  on  the  number  of  interfering  users  and  their  received 
power  at  the  base  station.  Therefore,  one  should  first  derive  the  bit  error 
probability  given  the  number  of  interfering  users;  from  that  one  should  derive 
the  packet  error  probability,  given  again  the  number  of  interfering  users,  and 
then  finally  average  the  latter  with  respect  to  the  Poisson  distributed  number 
of  interfering  users.  As  a  result,  the  best  that  can  be  done  using  this  method 
for  packets  (coded  or  uncoded  ones)  is  an  approximation.' 

5  NUMERICAL  RESULTS 

In  discussing  our  numerical  results,  we  must  note  several  things.  First,  we 
take  A  to  be  a  circular  region  of  radius  R  and  then  we  set  R  —  1.  This  can 
be  viewed  as  a  normalization  and  the  results  can  be  generalized  to  regions 
of  larger  or  smaller  size;  this  is  especially  evident  if  one  considers  that  the 
expected  number  of  users  in  A  is  given  by  N  =  A ttR2  and  that  our  results 
are  indexed  by  N  (and  therefore  by  A). 

A  brief  discussion  of  the  issue  of  normalizing  the  radius  R  of  the  cell 
is  in  order  here.  We  assume  that  in  the  mobile-to-base  mode  of  opera¬ 
tion,  all  mobiles  transmit  with  the  same  signal-to-noise  ratio  Eb/N o,  which 
is  such  that  the  signal-to-noise  ratio  received  at  the  base  from  a  mobile 
at  a  distance  r,  given  by  E'b/No  =  {Eb/No)  g(r),  is  at  least  as  large  as 
Eb/No  =  (Eb/No) g(R),  the  received  energy  from  a  mobile  at  the  periph¬ 
ery  of  the  cell.  Therefore,  if  Eb/No  is  set  at  the  minimum  acceptable  level 
for  reliable  communications,  the  transmitted  signal-to-noise  ratio  for  all  mo¬ 
biles  is  given  by  Eb/N0  =  (^Eb/No'j  / g(R );  if  R  =  1,  Eb/N0  =  "Eb/No, 
while  if  R  >  1,  g(R)  can  be  computed  and  used  to  compute  Eb/No  from  the 
above  formula,  while  using  R  =  1  for  all  other  calculations. 

It  is  also  important  to  consider  the  computational  requirements  of  our 
results.  Looking  at  (30)  or  (33),  one  can  see  that  there  are  multiple  levels 
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of  summations  and  integrals  involved  in  the  computation  of  each  value  of 
the  packet  error  probability  (note  that  Ef{ •}  is  an  integral).  In  computing 
the  results  shown  here,  we  sought  to  avoid  this  last  expectation  over  f  and 
explored  two  alternative  methods. 

The  first  method  was  to  set  f  equal  to  some  value  such  that  a  given 
probability  of  acceptable  communication  is  maintained.  As  discussed  in  [1], 
it  can  be  shown  that  in  this  case,  f  =  \/PR,  where  P  is  the  probability 
of  acceptable  communication  to  be  maintained.  If  one  considers  values  of  P 
close  to  one-as  would  be  the  range  of  values  one  might  like  to  examine-one 
can  see  that  the  results  obtained  will  be  overly  pessimistic. 

As  a  result  of  the  pessimistic  nature  of  the  first  method,  we  have  opted 
to  use  another  method,  namely  to  set  f  equal  to  a  value  such  that  an  average 
value  of  g(f)  is  obtained;  then  f  =  [22r{<7(r)}]~1/0',  where  a  and  <?(•)  are  as 
in  (2).  We  find  that  the  expectation  is  given  by 

Er  W)}  =  ^  +  (g  -  *2'")  =  S(f)-  (56) 

The  value  of  this  method  is  due  to  the  fact  that  the  probability  of  symbol 
or  packet  correctness  is  usually  a  convex  fl  function  of  g(r)  in  part  of  the 
range  of  r.  In  this  range,  we  can  use  Jensen’s  inequality 

-MPcGKr))}  <  Pc{Er(g(r))}  =  Pc(g(f))  (57) 

to  obtain  the  upper  bound  on  the  right-hand  side  which  depends  only  on 
g(r).  In  the  range  of  r  in  which  the  function  Pc  is  not  convex  ft  in  r,  we  can 
still  use  Pc(g(f)  as  an  approximation. 

Fig.  2  compares  the  results  using  the  Gaussian  approximation  tech¬ 
nique,  described  in  Section  III-A,  with  the  bound  presented  in  Section  III-B 
for  the  case  of  coherent  signaling  and  RS(32,  16)  coding  with  error-only  cor¬ 
rection  decoding.  This  demonstrates  that  the  bound  tends  to  be  extremely 
pessimistic.  On  the  other  hand,  the  Gaussian  approximation  tends  to  be 
moderately  optimistic  for  most  of  the  range  of  N  (the  average  number  of 
neighbors).  This  last  argument  is  justified  as  follows:  Gaussian  approxima¬ 
tions  tend  to  be  usually  optimistic  when  applied  to  multiple- access  problems 
(refer  to  [5]  and  [6]).  However,  the  assumption  that  all  interfering  users  cause 
full  hits  is  pessimistic.  Furthermore,  the  use  of  g(r )  results  in  pessimistic 
values  of  the  probability  of  correct  reception  for  part  of  the  range  of  r.  Thus 
overall  the  Gaussian  approximation  used  here  is  expected  to  be  moderately 
optimistic.  This  can  be  actually  verified  through  a  comparison  to  the  exact 
expressions  in  Fig.  6  and  7  below. 
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Fig.  3  deals  with  the  coherent  case  and  demonstrates  the  decrease  in 
throughput  with  decreasing  Reed- Solomon  code  rate.  As  one  might  expect, 
throughput  is  proportional  to  code  rate  due  to  the  normalization  by  the 
code  rate.  However,  the  RS(32,8)  code  provides  a  bit  more  than  half  of  the 
RS(32,16)  throughput  when  N  =  150,  which  is  explained  by  an  improve¬ 
ment  in  the  expected  probability  of  error.  In  both  Fig.  2  and  3,  the  values 
of  the  throughput  for  large  N  (larger  than  150)  are  optimistic  and  should 
not  be  considered  fully  dependable. 

We  now  move  on  to  the  noncoherent  case  (32-ary  FSK  with  noncoherent 
demodulation).  Fig.  4  corresponds  to  Fig.  2.  Once  again  it  is  demonstrated 
that  the  bound  is  very  pessimistic  and  the  Gaussian  approximation  overly 
optimistic  for  large  N . 

Fig.  5  presents  results  using  RS(32,  16)  and  RS(32,  8)  coding.  The  same 
observations  that  were  made  regarding  Fig.  3  should  be  made  here. 

Fig.  6  demonstrates  the  relationships  among  the  exact  expression  devel¬ 
oped  in  this  section  and  the  Gaussian  approximation  and  bound  presented 
in  Section  III.  The  bound  tends  to  be  very  pessimistic.  The  Gaussian  ap¬ 
proximation,  on  the  other  hand,  tends  to  be  somewhat  pessimistic  for  low  to 
moderate  values  of  N  (the  average  number  of  neighbors),  while  it  provides 
a  closer  approximation  and  even  becomes  slightly  optimistic  for  moderate 
to  large  values  of  N.  This  reflects  on  the  previous  throughput  results:  for 
small  N ,  the  approximation  is  closer  to  the  exact  value  than  the  bound  but 
is  pessimistic;  for  moderate  IV,  it  is  close  to  the  exact  value;  for  large  N , 
the  results  exceed  what  they  should  be  and  this  effect  is  magnified  by  the 
RS  coding. 

Fig.  7  shows  results  for  synchronous  and  asynchronous  systems.  These 
results  are  so  close,  that  the  appear  as  a  single  curve  in  the  figure.  The 
closeness  of  the  results  for  the  synchronous  and  asynchronous  cases  is  due 
to  the  fact  that  the  probability  of  a  hit  is  Ph  =  l/q  for  synchronous  systems 
and  Ph  =  (1  -f  1  /Nb)/q  for  asynchronous  systems,  whereas  m $  =  1/3  for 
asynchronous  systems  and  m $  =  1/2  for  synchronous  systems,  thus  balanc¬ 
ing  out  the  effect  of  other-user  interference  of  the  two  systems.  Also  shown 
is  a  curve  that  represents  an  asynchronous  system  with  Rayleigh  fading  (pa¬ 
rameter  b 2  —  0.01).  As  expected,  this  system  does  significantly  worse  than 
its  counterpart  without  fading. 
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6  CONCLUDING  REMARKS 

We  presented  methods  by  which  error  probability  and/or  throughput  can 
be  computed  for  FH  cellular  mobile  packet  redio  networks.  These  methods 
consist  of  a  Gaussian  approximation  applicable  to  all  systems  and  an  exact 
method  that  we  could  only  apply  to  coherent  systems.  It  was  shown  that 
the  previously  derived  lower  bound  on  system  throughput  (upper  bound 
on  error  probability)  is  a  loose,  pessimistic  bound  in  comparison  to  the 
approximation  and  the  exact  results.  The  approximation  was  shown  to 
have  satisfactory  accuracy  for  the  lower-to-moderate  values  of  the  average 
number  of  neighbors,  but  to  be  overly  optimistic  for  large  values  of  the 
average  number  of  neighbors. 

One  major  problem  encountered  lay  in  the  computational  complexity  of 
the  expressions.  In  general,  each  curve  in  the  figures  presented  required  a 
considerable  amount  of  CPU  time  to  compute.  In-the  case  of  some  of  the 
exact  expressions  (especially  the  asynchronous  case  without  fading),  this 
computation  was  extremely  long.  We  used  basically  a  VAX  780  computer 
and  a  few  hours  of  supercomputer  CPU  time.  What  this  indicates  is  that 
this  work  may  very  well  represent  a  practical  limit  as  to  how  precisely  one 
can  characterize  these  systems  and  still  be  able  to  compute  analytically  the 
various  performance  measures  of  interest. 
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Figure  1.  BFSK  SFH/SSMA  receiver  model. 


Figure  2.  Throughput  versus  Average  Number  of  Neighbors:  Gaussian  Approximation 
versus  Bound  (RS(32,16)  coding  with  error-only  correction  decoding,  coherent 
BPSK ,  SNR=1 5dB  ,  q  =  100,  N,=10) 


RS (32 , 16) 


Figure  3.  Throughput  versus  Average  Number  of  Neighbors  for  Varying  Code  Rate 
(Reed-Solomon  coding  with  error-only  correction  decoding,  coherent  BPSK, 
SNR=15dB  ,  q=l  00  ,  N,=10) 


Figure  A.  Throughput  versus  Average  Number  of  Neighbors:  Gaussian 
Approximation  versus  Bound  (RS(32,16)  coding  with  error-only  correction 
decoding,  noncoherent  32-ary  FSK,  SNR=15dB,  q=100,  N, =10) 


Figure  5.  Throughput  versus  Average  Number  of  Neighbors  for  Varying  Code  Rate 
(Reed-Solomon  coding  with  error-only  correction  decoding,  noncoherent  32-ary 
FSK,  SNR=1 5dB ,  q=100,  N,=10) 


Bound’ 


Figure  6.  Bit  Error  Probability  versus  Average  Number  of  Neighbors:  Exact 
Method,  Gaussian  Approximation,  and  Bound  (coherent  BPSK,  asynchronous, 
SNR=1 OdB ,  q=l 00 ,  N,=10) 


